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ABSTRACT 

The Swift satellite early x-ray data shows a very steep decay in most of the Gamma-Ray 
Bursts light curves. This decay is either produced by the rapidly declining continuation of 
the central engine activity or by some left-over radiation starting right after the central engine 
shuts off. The latter scenario consists of the emission from an "ember" that cools via adiabatic 
expansion and, if the jet angle is larger than the inverse of the source Lorentz factor, the large 
angle emission. In this work, we calculate the temporal and spectral properties of the emission 
from such a cooling ember, providing a new treatment for the micro-physics of the adiabatic 
expansion. We use the adiabatic invariance of p\ /B (p± is the component of the electrons' 
momentum normal to the magnetic field, B) to calculate the electrons' Lorentz factor during 
the adiabatic expansion; the electron momentum becomes more and more aligned with the 
local magnetic field as the expansion develops. We compare the theoretical expectations of 
the adiabatic expansion (and the large angle emission) with the current observations of the 
early x-ray data and find that only ~ 20% of our sample of 107 bursts is potentially consistent 
with this model. This leads us to believe that, for most bursts, the central engine does not 
turn off completely during the steep decay of the x-ray light curve; therefore, this phase is 
produced by the continued rapidly declining activity of the central engine. 
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1 INTRODUCTION 

The central engine of the GRBs (Gamma-Ray Bursts - see Piran 
2005 for a review) is hidden to direct observations and its workings 
are largely unknown. The only information that we currently have 
about the GRB is obtained from its electromagnetic radiation. We 
have to look for signatures in the radiation mechanism to under- 
stand how Nature produces these outbursts. 

Swift has provided very early x-ray data that shows that for 
most bursts there is a very steep decay lasting for about 10 minutes 
(Tagliaferri et al. 2005, Nousek et al. 2006 - see Zhang et al. 2006 
and references therein for possible physical explanations). These 
observations suggest that the rapidly declining x-ray light curve 
(LC) and the burst are produced by the same source, because the x- 
ray LC, when extrapolated backwards in time, matches the gamma- 
ray LC (O'Brien et al. 2006). Therefore, a natural question arises: 
Is the early x-ray data really just a rapidly declining continuation of 
the central engine activity, originally seen in the gamma-ray band, 
but now seen at lower energy? Or does the central engine switch 
off abruptly and the early x-ray data doesn't reflect the activity of 
the central engine? 
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If the central engine completely shuts off when the gamma- 
ray photons' flux falls below the gamma-ray detector sensitivity, 
then the emission from a cooling "ember" would be responsible for 
the early x-ray steep decay. This source - which had just produced 
the gamma-ray emission - would be cooling by adiabatic expansion 
(AE). In this paper, we study the flux properties of a "hot" shell that 
undergoes AE and cools. The goal is to determine if the observed 
x-ray steep decay data is consistent with the AE scenario. If so, 
then the central engine did shut off abruptly right after the gamma- 
ray emission ceased. On the other hand, if it isn't, then the data is 
reflecting the rapidly declining activity of the central engine. The 
reason for this is that any other process that does not invoke central 
engine activity to explain the early x-ray data has problems explain- 
ing the smooth temporal connection observed in the LC between 
the prompt emission and the early x-ray data. 

AE has been studied before to predict the long-wavelength af- 
terglow from GRBs (Meszaros & Rees 1997) and also, the optical 
flashes from internal and reverse shocks (Meszaros & Rees 1999, 
Sari & Piran 1999a). In this work, we describe the evolution of a 
collisionless plasma due to AE and show that this is in general dif- 
ferent from AE of an ideal gas. 

We first present the micro-physics of the AE for a collision- 
less plasma in §2. Then, we use it to calculate the flux properties of 
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a source undergoing AE in §3. We look at what the current obser- 
vations tell us in §4, and we put them in the context of the central 
engine in §5. We summarize our results and give our conclusions 
in §6. 



2 MICRO-PHYSICS OF THE ADIABATIC EXPANSION 

For an ideal gas, the pressure evolves due to adiabatic expansion as 
P e j oc p a " oc V~ a " , where a e = 4/3 for a relativistic gas and p 
and V are the co-moving density and ejecta volume, respectively. 
The collisions between electrons are extremely rare in GRB rela- 
tivistic shocks, therefore, one needs to be careful in the use of this 
formul^B 

For a collisionless magnetized plasma, assuming that no col- 
lective plasma processes randomize the particles' velocity due to 
scattering, particles move along a magnetic field line and, by using 
the concept of adiabatic invariant (Jackson 1998, Rybicki & Light- 
man 1979), we can calculate the particles' energy. This invariant 
describes that, for slowly varying fields, the magnetic flux through 
the orbit of the particle is a constant, or p\/B is an adiabatic in- 
variant, where p± is the component of the particle's momentum 
transverse to B, the magnetic field. For highly relativistic particles, 
p± ~ m e c7x, so that 

Jx/B = constant (1) 

can be used, where 7x is the Lorentz Factor (LF) of the electron 
in the transverse direction, m e is the electron's mass and c is the 
speed of light (from now on, c = 1). This relationship can be used 
because the magnetic field decays on a much larger length-scale 
than the electron's gyro-radius (see Appendix A). It is worth noting 
that the parallel component of the electron's momentum remains 
unchanged; this will be briefly discussed in the last section. 

In the next sections, we will make use of QJ to predict the 
evolution of the electrons' LF in a hot shell that undergoes AE. 
We will use it to calculate the properties of its synchrotron and 
synchrotron-self-Compton (SSC) radiation. 



3 ANALYTICAL LIGHT CURVES OF AN 
ADIABATICALLY COOLING EMBER 

Let's assume that the GRB ejecta was heated by some process 
(shocks or magnetic dissipation) and suddenly the central engine 
switches off completely. There is no other energy injection mecha- 
nism at hand, so it begins to coast (the LF of the ejecta is constant, 
see §3.6) and cools via AE. We will calculate the flux properties of 
this cooling ember. 

1 We have estimated the mean free path for Coulomb scattering between 
a hot electron and a cold electron and find it to be much larger compared 
to the shell thickness. The electric field associated with a relativistic hot 
electron is calculated using the Lienard-Wiechert potential, and we find the 
cross-section for a significant interaction, i.e. leading to a fraction of the en- 
ergy of the hot electron transfered to a cold electron, is close to the Thomson 
cross section (a ss or/3). The mean free path is A = (ntr) , where n 
is the electron density. From the total energy E, the distance of the shell 
from the center of explosion R, the source LF T, and the co-moving shell 
width A', we find n = E / (Tm p c 2 4:TrR 2 A'). Using the usual notation 
Q n = Q/10", we obtain: A/A' ss 700-Rf 5 r 2 -E~ 2 1 . For scattering be- 
tween hot electrons, the conclusion is the same. Therefore, Coulomb scat- 
tering between the electrons is not significant. 



3.1 Ejecta width and magnetic field 

In the following subsections we will provide the basic ingredients 
for the radiation calculation. First, we need to determine the co- 
moving thickness of the ejecta, which could be: 

A' = R/F or A' = A . 

We will call these cases: thin ejecta shell (an ejecta that undergoes 
significant spreading) and thick ejecta shell (an ejecta that experi- 
ences no significant spreading), respectively. The observed time is 
t oc R/T 2 , where 7? is the distance of the source from the center of 
the explosion and V is the LF of the source with respect to the rest 
frame of the GRB host galaxy. 

The magnetic field in the GRB ejecta could be a combination 
of frozen-in field from the central explosion and field generated lo- 
cally. We prescribe the decay of the field by using the flux-freezing 
condition (Panaitescu & Kumar 2004), which gives: 

B± oc {RA'y 1 and B n oc R~ 2 . 

We will use the field that decays slower, although this is highly 
uncertain. This is because the magnetic field generation mechanism 
for GRBs is still not well understood, so the relative strength of B± 
and 73 1 1 is unknown. 

3.2 Electrons' energy distribution 

For an adiabatically expanding source, no more energetic electrons 
are injected in the system when the shock has run its course. This 
means that, after some time, there will be few electrons with ener- 
gies higher than the cooling electron LF, j c . Therefore, the electron 
population above 7 C will be truncated due to radiation losses and 
the emission for v c < v will rapidly shut off (v c is the synchrotron 
frequency corresponding to electrons with 7 C ). At this point, the 
electron distribution will follow oc y~ p for 7* < 7 < 7 C , where 
ji is the typical LF of the electrons, since we would be dealing 
only with adiabatic electrons. Moreover, since the radiative cooling 
quickly becomes less important than the adiabatic cooling because 
the magnetic field decays rapidly with the expansion of the ejecta, 
both 7; and y c will evolve in the same way. 

For the case -y c < 7,, the radiation losses would dominate and, 
after some time, they would make the whole electron distribution 
collapse to a value close to j c . A narrower range in the electron 
distribution would be responsible for the radiation. In this paper we 
focus on the 7, < 7 C case. 

3.3 Basics of Synchrotron and SSC 

The electrons' four-momentum is given by P = 
m e (7, 7 cos «', 7 sina', 0), which can be also written as 
P = m e (7,7||,7x,0), where 711 is the component of the elec- 
trons' momentum parallel to B. The pitch angle, which is the 
angle between B and the velocity of the electrons, is a' . In this 
notation, the electron's LF is then 7 2 = 7y + 7x- According to 
the prescription of the adiabatic invariance, 7x evolves following 
(1) and 7|| remains unchanged. We assume that at the onset of the 
adiabatic expansion 7x ~ 711, then quickly when time doubles, 
the radius would have also doubled, making the magnetic field 
decrease by at least that factor and reducing 7x making 711 > 7x, 
which gives: 



7 = 7|u/l + Hf ~ Til and Sina = — ~ — . (2) 

V 7 il 7 7 " 
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As the transverse component of the momentum decreases, due to 
the decay of the magnetic field, the pitch angle decreases, which 
makes the electron's momentum more aligned with the local mag- 
netic field. 

Knowing the electrons' energy distribution, the emission at 
any given frequency and time can be calculated using the syn- 
chrotron spectrum: 



{v/v a ) 2 (v a /vif^ v<v a 

(v/Vi) 1/3 Va <V <Vi 

(iz/i/i)-^- 1 ^ 2 Vi<V< v c , 



(3) 



for the case v a < v t < v c , where v a is the self absorption fre- 
quency and it is obtained using equation (52) of Panaitescu & Ku- 
mar (2000) (see, e.g., Katz & Piran 1997, Sari & Piran 1999b). The 
characteristic synchrotron frequencies are obtained from the corre- 
sponding electrons LFs: 



(4) 



27rm e (l + z) 

where 7i, c and sin a' are given by (2). The observed peak flux is 

„ (1 + z)V3e 3 N e BF . , 

~ ±nd 2 (z)m e SmQ ' (5) 

where (II is the luminosity distance, N e is the number of radiating 
electrons (which in this scenario is constant), z is the redshift and 
e is the electron's charge. 

For the SSC case, the flux peaks at v^ 2 with magnitude 
r e F Vi , where r e = N e ar /(4irR 2 ) is the optical depth to elec- 
tron scattering. We will calculate SSC emission for photons above 
v a . We will use the same synchrotron piece- wise spectrum, which 
is just a very crude approximation. 



3.4 Temporal and spectral properties 

For synchrotron emission of a cooling ember undergoing AE, we 
find: 



F Vi oc t- 3 (t- 3/2 ), 
and 



ott- :i (t- :i/2 ), Va ott- L2,b {t-» ,b ){6) 



12/5,.-9/5n 



t 2 (t 2 y 

r v < t-*(t->i/3 

t -3(p+l)/2^-3(p+l)/4^-(p-l)/2 



V < V a 

V a < V < Vi 

Vi < V < v c , 



(V) 



where the time dependences are reported for the thin ejecta shell 
and parenthesis are used for the thick ejecta shell. 

Using the same notation as above, for SSC emission, we find: 



and 

F v oc 



t -(3 P +7)/2 (r (3 P +ll)/4 )i/ -( P -l)/2 Vil l<y<y cl l 



(9) 



One can see that for the synchrotron case, the flux decays 
rapidly for Vi < v < v c ; for the SSC case, the flux decays even 
more rapidly for Vij 2 <u< v c ^ 2 . For both cases, the peak frequen- 
cies of the spectrum also show a fast decrease with time. Also, for 
both cases, the thin ejecta case gives a faster time decay, since the 
shell spreads significantly, allowing the ejecta to cool faster. 

To compare our theory with the observations, we provide re- 
lations between the temporal decay index a and the spectral index 
(3 in Table 1, using the convention F v oc i~ c V~' 3 . 





Synchrotron 

Vi < v < v c 


SSC 

va 2 <v< v c j 2 


Thick ejecta 






a = 1.5/3 + 1.5 


a = 1.5/3 + 3.5 


(A' = A ) 






Thin ejecta 
(A' = R/T) 


a = 3/3 + 3 


a = 3/3 + 5 



Table 1. Closure relations between a (decay index) and /3 (spectral index) 
for a cooling ember undergoing AE (to = t c ). 



To summarize, the emission from an adiabatically cooling 
source has the following properties: 

(i) Its spectral index must be equal to the one at the end of the 
prompt emission phase of the gamma-ray burst, /3 7 . 

(ii) The temporal decay index must obey one of the closure re- 
lations in Table 1 . 

(iii) The peak frequency of the spectrum should decrease with 
time as predicted in (6) and (8). 

(iv) After some time, on the order of t c (defined in §3.5), the 
spectrum should have an exponential cut-off at frequencies greater 
than the cooling frequency (§3.2). 

Points (ii) and (iii) have to correspond to the same radiation 
mechanism (synchrotron or SSC) and the same ejecta width case 
(thin or thick). 

If one were to consider the electrons' energy as given by the 
adiabatic expansion of an ideal relativistic gas, instead of using the 
methods of adiabatic invariance, then 7 oc V~ 1/3 (§2, see also §3 
of Meszaros & Rees 1999). In this case, the velocity distribution 
of the electrons during the adiabatic expansion phase is isotropic, 
therefore, sin a' is a time independent constant of order unity. We 
calculate the temporal decay indices as done above and find the 
following results. For synchrotron: a = 2.3(3 + 1 and a = 4/3 + 2, 
and for SSC: a = 3.7(3 + 3 and a = 6(3 + 4, for the thick and 
thin ejecta cases, respectively. The difference in the temporal decay 
indices for the synchrotron case compared to the ones on Table 1 is 
< 20% for (3 e [0.5 — 2]. Because SSC has a stronger dependence 
on 7, the difference we find in a is larger. 

So far, we have considered only the flux-freezing condition 
to prescribe the evolution of the magnetic field, but we can also 
determine the magnetic field using the equipartition consideration, 
i.e. the energy density in the magnetic field is a constant fraction 
of the electrons' internal energy density (Sari, Piran & Narayan 
1998). To obtain this last quantity, one needs to know 7, which 
could be obtained either using the adiabatic invariance methods or 
by the ideal gas law - as mentioned in the last paragraph. But for 
the equipartition consideration we will only consider the ideal gas 
law, because if there is a mechanism that maintains an equiparti- 
tion between magnetic energy and electron energy, then that same 
process is also likely to keep different components of electron mo- 
mentum coupled and that will lead to an ideal gas expansion law 
for electrons. Therefore, the magnetic field in this case is given by 
B 2 oc V~ 4/3 , where V oc R 2 A'. The synchrotron and SSC emis- 
sion decays, for the thick case, are both steeper by 0.3(3 + 0.3 than 
the ones presented on the last paragraph, but both thin cases remain 
unchanged. 
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3.5 Large Angle Emission 

If the central engine switches off abruptly and the gamma-ray pro- 
ducing ejecta has a opening angle 9j, such as 9j > Y^ 1 , Large 
Angle Emission (LAE) will be also present (Fenimore & Sum- 
ner 1997, Kumar & Panaitescu 2000). The LAE flux declines as 
a = 2 + (3 and the peak frequency of the spectrum decays as t~ . 
Therefore, the AE flux generalljQ decays faster than LAE's and the 
AE's peak frequencies always decrease faster than LAE's. 

The time-scales for these two phenomena, LAE and AE, are 
essentially the same, they are set by 

*c=2^. (10) 

LAE and AE start at the same time, to , and same site, Ro : right after 
the central engine has switched off, and the fluxes decline with time 
as: 



(11) 



where Fo is the flux at to and a is the decay index of either LAE 
or AE. The shape of the LC depends on the values for to and t c 
(Figure 1); for §3.4, to = t c . The case to < t c is unphysical, since 
it implies that when AE starts, the shell's electrons have not yet 
cooled substantially, i.e. the shell's radius hasn't doubled. 

If Oj < r _1 , then there will be no LAE, so AE would be the 
only emission present after the central engine turns off. On the other 
hand, if 9j > T~ , then LAE will dominate over AE (see footnote 
2). LAE will cease with the detection of the last photons coming 
from 6j and, at this time, the flux will smoothly become dominated 
by the AE emission, i.e. there will be a break in the LC to the power 
law decay for AE (Figure 1: Bottom). The photons from 9j will 
arrive at a time tj « to + R9 2 /2 =t + 9 2 Y 2 t c . 



3.6 Electron-positron pair-enriched ejecta 

When the ejecta cools by adiabatic expansion, the thermal energy of 
the protons and electrons is converted back to bulk kinetic energy 
of the shell, increasing F. Even in an extreme case where all the 
electrons' energy goes into the shell expansion, V increases onl 
by a factor of ~ 2, if the protons and electrons energy is < m p c' 
Therefore, the effect of this change to the observed flux is less than 
a factor of 2, a relatively small effect. For this reason, we have used 
a constant Y for the calculations done so far. 

On the other hand, if the ejecta consists of e pairs, then the 
increase in Y during the adiabatic expansion would be considerable, 
and it would scale as oc j^ 1 (the observed energy in the shell is a 
constant and scales as oc 7F). Since F increases, then the observed 



time is t — to = J } 



F v = Foil — 



Ro 
to 



3t c 



dR/(2Y ). For the thin ejecta case, we find: 



(12) 
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Figure 1. The normalized flux density, equations (11) and (12), plotted vs. 
observed time, assuming that the observed frequency v is always between 
Vi and v c (if u c < u, then there is no AE, only LAE if 0j > see 
§3.2). This emission is produced by the last ejected shell, because con- 
tributions from previously ejected shells would be buried in the emission 
of subsequent shells, since both LAE and AE decay very fast. Top. Using 
to = 100s and t c = 10s. The LAE and AE-Baryonic decay indices corre- 
spond to a = (3,6), respectively, and the AE-Pair has S = 8 (/3 = 1, using 
AE: Synchrotron - thin ejecta). Bottom. Using to = t c = 100s and the 
same a's and S as above. In this illustrative example we have Oj = 2/T, 
therefore, LAE dominates over AE until tj = 500s, when AE-Baryonic 
takes over. This break in the LC from the LAE to the AE power law decay 
(which should be a smooth transition and it is done in the figure for display 
purposes) has never been observed. 



2 Except when /3 < 1 for the AE case of synchrotron emission from a 
thick ejecta shell. 

3 We have assumed a co-moving observer sitting in the middle of an infinite 
parallel shell that sees the left and right halves of the shell move away from 
him. Assuming the electrons' LF in the shell rest frame is 10 3 (and that 
the protons are essentially cold since the heating mechanism energized all 
particles equally), the LF of the shells would be 1 + 10 3 m e /m p = 1.5. 
An observer far away would mainly detect radiation from the half moving 
towards him, since the radiation is beamed. This observer would see that 
the LF of this half has increased by a factor of ~ 2. 



with characteristic frequencies oc [1 — (t — to)/(3t c )] 4 , where 
5 = (4f3 + 4, 4/3 + 10) for the synchrotron (ui < u < v c ) and 
SSC (vi'y'i VcJ 2 ) cases, respectively. The time decay index, 
for the to = t c case, is 8 puffin ■ therefore! the LC steepens con- 
tinuously. It decays even faster than the AE-Baryonic case, since 
the observed time gets compressed because Y is increasing consid- 
erably. If 9j > r~\ then LAE prevents it from steepening more 
than 2 + 13: completely taking over the emission since essentially 
to (Figure 1: Bottom). 
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3.7 Reverse Shock emission 

In this short subsection, we explore the possibility that the GRB 
ejecta, that just produced the prompt emission, interacts with the 
interstellar medium (ISM) and a reverse shock (RS) crosses it. The 
ejecta cools adiabatically after the RS has passed through it and we 
assume that it follows the Blandford-McKee self similar solution 
(Blandford & McKee 1976), during which V decays in time. 

Using the same methods as in the previous subsections, we 
can calculate the LF of the electrons in the ejecta after the passage 
of the RS (see Appendix B). We find that the RS flux decays as 
oc t- ill(p+1)/56S = t~ 2 53 (synchrotron emission: u t < v < v c , 
for a thin shell, using B± and p = 2.5), which gives the clo- 
sure relation: a = 1.45/3 + 1.45 (and a = 1.45/3 + 1.67 for 
SSC). If we determine the electrons' LF using the ideal gas law, 
then the RS synchrotron flux would decay as oc £-( 20 p+ 7 )/ 24 = 
t -238 (for the same case as above), which is still steeper than the 
oc f-( 7 3p+2i)/96 = £-2.12 derived by Sari & p iran (1999^ w h e re 

they used 7 oc V -1 ' 3 and the equipartition consideration. 



4 APPLICATION TO THE GRB EARLY "AFTERGLOW" 

AE, together with LAE, dictates the emission of the source after the 
central engine has completely turned off. In this section, we will 
determine if the early x-ray steep decay observed by Swift obeys 
our theoretical LCs for AE and LAE0. We will do this for each one 
of the cases considered in the previous section. 

The early x-ray data shows a single power law decay with 

3 < a < 5 (Nousek et al. 2006, O'Brien et al. 2006, Willingale et 
al. 2007). With this information, the to > t c case can be ruled out, 
since, for this case, the theoretical shape of the LCs for LAE and 
AE-Baryonic is inconsistent with the early x-ray observations and 
the AE-Pair LC decays extremely fast (Figure 1: Top). Therefore, 
we focus on the to = t c case only. 

The next possibility we explore is to see if the early x-ray data 
obeys LAE or AE (from a baryonic ejecta). To check the validity 
of these two scenarios, respectively, we will take a sample of bursts 
and see how many cases are possibly consistent with either LAE or 
AE. 

Our sample consists of 107 GRBs for which their spectral in- 
dex during the early x-ray decay (J3 X ) and their temporal decay in- 
dex during this phase have been previously determined (the sam- 
ple of Willingale et al. 2007). We first select the bursts for which 
/3 7 = (3 X , which cuts down the sample to 55 bursts. Eight of these 
bursts show strong spectral evolution, inconsistent with LAE and 
with AE (Zhang et al. 2007: Zhang et al's sample essentially con- 
tains all our sample), which leaves us with 47 bursts. Moreover, 
we check how many of these satisfy the relations between a and 
/3 for LAE or AE (Table 1) within about a 90% confidence level, 
and that narrows down the sample to 20 bursts. In conclusion, only 
a small percentage of the sample, 19%, is consistent with LAE or 
AE, which leads us to suggest that, for most bursts, the early x-ray 
data results from some other process, and the most natural conclu- 
sion is continued central engine activity. 

4 If the prompt emission is attributed to synchrotron, then, for some frac- 
tion of the parameter space, the radiative cooling timescale, t ra£ j is less than 
t c . However, shortly after the onset of the adiabatic expansion, t rac [ > t c , 
since the magnetic field decays rapidly with the expansion of the ejecta. For 
the SSC case, t ra d > t c is very likely at the onset of the adiabatic expan- 
sion, making the radiative cooling unimportant. 



It has also been claimed that the 7-ray emission extrapolated 
to x-ray energies, together with the early x-ray data, can be well 
fitted with a falling exponential followed by a power law (O'Brien 
et al. 2006). At first, this could be thought to be explained by a pair 
ejecta with 9j <r _1 undergoing AE, since its LC also steepens 
continuously (Figure 1: Bottom). However, for 3 bursts that show 
this continuous steepening: GRB 050315 (Vaughan et al. 2006, 
Lazzati & Begelman 2006), GRB 050724 (Barthelmy et al. 2005) 
and GRB 060614 (Mangano et al. 2007), the theoretical LC decays 
too fast and can't fit the observed early x-ray LC. Therefore, we 
rule out the possibility that the observed early x-ray decay is from 
a pair ejecta undergoing AE. 

Finally, we use our sample to test if the observed early x-ray 
steep decay is consistent with the closure relations derived for the 
RS (§3.7). In this case, the condition /3 7 = (3% is not necessary, 
therefore we start with the entire sample: 107 bursts. We eliminate 
19 of these, which show strong spectral evolution, inconsistent with 
RS. Out of the remainder, only 26 bursts (24%) are possibly con- 
sistent with the RS (16 of these are simultaneously consistent with 
LAE and AE). However, this scenario could be ruled out since it 
would be difficult to explain the connection observed in the LC 
between the prompt emission and the early x-ray data. The only 
way they could be smoothly connected is if the prompt gamma-ray 
emission is produced also by the RS, for which we lack evidence. 



5 DISCUSSION: THE CENTRAL ENGINE 

The results of the last section lead us to believe that the observed 
early x-ray decay for > 70% of GRBs is produced by the rapidly 
declining continued activity of the central engine, if we assume that 
there is a one-to-one correspondence between the temporal behav- 
ior of the central engine activity and the observed emission (Figure 
2). 

There are two other arguments that support the idea that, for 
most bursts, neither LAE nor AE might be consistent with the ob- 
served early x-ray decay. First, for some bursts, a break frequency 
has been seen passing through the x-ray band during the early steep 
decay, and it evolves faster than the oc t _1 expected in LAE: oc t~ 2 
for GRB 060614 (Mangano et al. 200"fl) and oc t~ 4 ~ t~ 3 for 
GRB 060904A (Yonetoku et al. 2008). Second, if AE is entirely 
responsible for this phase, then 9j has to be very small. If this is the 
case, then we should have observed the edge of the jet (a jet break 
with a oc t~ p optical LC) very early on. After inspecting many 
optical LCs (Butler & Kocevski 2007, Liang, E.-W. et al. 2007, 
Liang, E.-W. et al. 2008, Melandri et al. 2008, Panaitescu & Ves- 
trand 2008), we can conclude that most bursts with available early 
optical data don't show this expected jet break starting at very early 
times, i.e. t < a few hours. 

One way that we might have missed these early jet breaks 
could be explained by the "porcupine" model. In this model, the 
central engine ejects many very small angle (8j < T _1 ) jets. One of 
these jets, directed towards the observer, produces the gamma-ray 
emission. The central engine shuts off, and then the AE emission 
is produced. After a short time, all these small jets combine and 
give rise to a single jet with 9j > that interacts with the inter- 
stellar medium, giving rise to a forward shock (optical afterglow). 



5 Mangano et al. mention in their work that this could be attributed to the 
v c decrease due to adiabatic cooling of shock heated shells after an internal 
shock. 
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This model leaves no sign of a jet break. Another possible expla- 
nation for the lack of a very early jet break is that there is energy 
injection to the forward shock, making the oc t~ v optical LC more 
shallow. This last scenario is unlikely, as a large amount of energy 
is required - more than a factor of 10 increase - to make a oc t~ p 
LC as shallow as ~ t , which is the usual observed optical LC 
decay. 

It has been suggested that the observed early x-ray decay is 
produced by the forward shock (FS) driven by the ejecta interact- 
ing with the ISM (Panaitescu 2007). This scenario has problems 
explaining the smooth temporal connection in the LC between the 
prompt emission and the early x-ray steep decay. 

The model presented in this paper has several uncertainties. 
For example, we are not specifying how the magnetic field is gener- 
ated during the prompt emission. We are assuming: (i) that the mag- 
netic field coherence length-scale is larger than the electron gyro- 
radius, based on the observations (see Appendix A), and (ii) that 
there are no collective plasma effects that randomize the electrons' 
velocity. These assumptions allow us to use the adiabatic invariant 
presented in §2. If, for some reason these conditions are violated, 
then the adiabatic invariance of electron magnetic moment can't be 
used. However, the light-curve from an expanding shell is similar 
whether we follow electron cooling via adiabatic invariance or ideal 
gas law and therefore, the main conclusions we have presented here 
are unchanged. Another assumption made in the application of this 
model is that Swift's x-ray telescope band lies between Vi and v c (in 
that order), which can be inferred from the spectral information of 
most bursts during this phase. If, however, v a < Ui, then there will 
be no emission coming from the part of the shell that lies within an 
angle of T - to the observer line of sight - the only emission would 
be from LAE (if 0j > L _1 ). 



6 CONCLUSIONS 

We have explored the situation in which the central engine shuts off 
and the ejecta cools via adiabatic expansion (AE). We have derived 
and discussed this emission's temporal and spectral properties us- 
ing a new treatment for the micro-physics of the AE: the adiabatic 
invariant J±/B, describing the electron momentum normal to the 
magnetic field (see summary in §3.4). At the onset of the adiabatic 
expansion, as B decays rapidly, this component of the momentum 
decreases while the parallel one remains unchanged, making the 
electrons more and more aligned with the local magnetic field as 
the ejecta expands. The adiabatic invariant enables us to calculate 
the electrons' energy for a collisionless magnetized plasma, if no 
other collective plasma effects that randomize the electrons' veloc- 
ity are present. 

In regards to the central engine activity, we can draw a con- 
clusion: The fastest way that the observed flux can decline after the 
central engine shuts off is set by the Large Angle Emission (LAE) 
and the Adiabatic Expansion cooling (depending on the value of 

Oil 

The early x-ray steep decay shown in most of Swift bursts has 
been attributed to LAE. In this paper, we consider both AE and 
LAE for the very early x-ray data. LAE and AE both start with 
the assumption that the central engine shuts off abruptly. Only ~ 
20% of our sample of 107 bursts is possibly consistent with either 
LAE or AE, thereby suggesting that the observed early x-ray steep 
decay for a large fraction of bursts might be produced by the rapidly 
declining continuation of the central engine activity. 

The component of the electron's momentum parallel to the 
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Time 



Figure 2. Two possible scenarios for the contribution of LAE (dotted), AE 
(dashed) and the rapidly declining central engine activity (solid), as seen in 
the x-ray band. The two "humps" represent the 7-ray detection (from Swift 
BAT), extrapolated to the x-ray band, attributed to activity of the central 
engine. Top. The case for which the central engine activity drops extremely 
fast and LAE and AE appear. Only ~ 20% of our sample is possibly con- 
sistent with this scenario. Bottom. Our preferred scenario, where the central 
engine activity is the dominant contribution and it decays slower than the 
theoretical LCs for LAE and AE. 



magnetic field is unconstrained by the adiabatic invariance. This 
component would probably cool via Inverse Compton scattering 
with synchrotron photons. Another possibility is that the electrons 
are scattered by small scale fluctuations in the magnetic field, which 
would effectively couple the parallel and perpendicular compo- 
nents of their momentum, resulting in an adiabatic cooling similar 
to that in the ideal gas case, §2 (personal communication, Granot). 

Any process that does not rely on the central engine activity to 
explain the observed early x-ray steep decay (i.e. RS, FS) has prob- 
lems explaining the smooth temporal connection observed in the 
LC between the prompt emission and the early x-ray steep decay. 
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APPENDIX A: ELECTRON GYRO-RADIUS VS. 
MAGNETIC FIELD LENGTH-SCALE 

The observed peak energy in the prompt phase of a GRB is given 
by (4) and it is u, = (1.15 x 10~ 8 eV)B7 2 r(l + z)" 1 , assuming 
synchrotron emission. The value of V can be constrained to be a 
few hundred. There is a wide range of allowed values for 7 i. For 
7* = 10 3 , Vi = lOOfceV, T = 100 and 2 = 1, then B = 2 x 10 5 G, 
and in that case, the electrons' gyro-radius is r — m e c 2 ji/(eB) w 
10cm. 

If the magnetic field responsible for the prompt emission is 
the frozen-in field from the central explosion, then it decays on a 
length-scale on the order of the source size, which is much larger 
than the electrons' gyro-radius and the adiabatic invariant presented 
in § 1 can be used. 

If the field is produced locally (e.g. by the Weibel instability), 
then we need to estimate its coherence length and compare it with 
the electrons' gyro-radius. For instance, magnetic field generated 



by the Weibel instability will have a coherence length on the order 
of the plasma length As = c/uj p , where uj p = (47re 2 n/m e ) 1 ^ 2 = 
6 x W i (n) 1/2 s' 1 and n is the co-movina electron number den- 
sity in units of cm -3 . For the prompt emission, recent studies have 
shown that the radius of emission is on the order of 10 15-16 cm or 
even larger (Kumar et al. 2007, Racusin et al. 2008, Zou, Piran & 
Sari 2009, Kumar & Narayan 2008), therefore n = 5 x 10 5 ~ 8 cm _3 
(see footnote 1), which gives As = 20 — 700cm, making As larger 
than r by at least a factor of 2. Moreover, this magnetic field decays 
extremely fast in time (in about ujp 1 ), which would be less than 
10 _s s in the source co-moving frame or 10~ 10 s in the observer 
frame. This locally generated magnetic field cannot be responsible 
for the prompt emission, unless it is sustained for at least ~ Is (the 
co-moving time-scale of a few mili-second prompt pulse), which 
would require a much larger As (Keshet et al. 2008 mention that 
the field must persist over 10 10 As downstream, which in this case 
would be ~ 10 12 cm). Therefore, it is safe to assume that even if the 
field is generated locally r « As, allowing us to use the adiabatic 
invariant. 

The prompt emission phase could also be attributed to the SSC 
emission, which requires smaller values for 7 j and B. For 7 j = 
10 2 , then B = 2 x 10 3 G, which gives r « 100cm ~ As. But 
as mentioned above, the field has to be coherent on length-scale 
~ 10 12 cm in order that it does not decay away on time < Is. 



APPENDIX B: REVERSE SHOCK EMISSION 
CALCULATION 

After the RS has crossed the ejecta, it follows roughly the 
Blandford-McKee self-similar solution (Blandford & McKee 
1976), in which the bulk LF and pressure of the shocked ISM are 
given by: 

7 (t,r) =l(t) X ~K P(r,t) =4m p c 2 n[ 7 (t)] 2 x"^, (B-l) 

where y(t) oc J£ - ( 3-s )/ 2 is the LF of material just behind the 
shock, n oc R~ s is the ISM particle number density, \ is the simi- 
larity variable and m p is the proton mass (see, e.g, Sari 1997). 

Let us assume that the ejecta is at Xej and it has a pressure P ej 
and a LF V e j, which - because of pressure and velocity equilibrium 
at the contact discontinuity - should be the same as the bulk LF and 
pressure of the shocked ISM at Xej (B-l). The pressure in the ejecta 
for the thin case (A' = R/T e j) is given by P e j oc V f±, where 

V oc R 2 A', and using the adiabatic invariance (1), we can use 

1/2 

7_l oc _B X , where B± is given in §3.1. Therefore, we have P e j oc 
7T 4 r 3 f and using (B-l) and the contact discontinuity equilibrium 
conditions we obtain: 

We can solve for r ej in terms of R, which gives V e j oc 

2(63-32s + 4s 2 ) 

32-73 _ p or uniform ISM (s = 0), we can determine 
the observed time by t — to = J R dR/(2T 2 j) and following the 

procedure on §3.3, we can find that the RS shock flux decays as 

oc t-*Kp+i)/568 = t -a.63 for p = 25 (u . <v< Vc) 

We can also determine the electrons' LF by using the ideal 
gas law. For this case, 7 oc V _1,/3 , therefore the pressure in 
the ejecta is P e j oc V _4//3 . For the thin ejecta case, we can re- 
peat the calculation done before (B-2), just modifying P e j, and 
we obtain V e j oc R~^ 7 ~ 2a ^ 2 . To calculate the RS synchrotron 
emission for uniform ISM, we can use B± oc (i?A') _1 and we 
can follow the procedure on §3.3, but as mentioned before, for 
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this case sin a' ~ 1. We find that the RS decays in this case as 
oc £-(20p+7)/24 = t -2.38 for p = 2 5 (y. <v< Vc y if we deter- 
mine the magnetic field using the equipartition consideration and 
determine the electrons' LF with the ideal gas law, then we get the 
same result reported by Sari & Piran (1999a). 



